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SUMMARY 


•He  well  known  chemical  equilibrium  problem  la  expresaed 
In  the  form  of  minimizing  the  free  energy  of  a  mixture  In 
order  to  compute  the  chemical  composition  at  equilibrium. 

By  piece-wise  linear  approximations  to  the  free  energy 
function,  the  problem  becomes  a  linear  program  which  can  be 
solved  by  a  standard  code  on  a  computing  machine.  Successive 
approximations  give  any  degree  of  accuracy. 
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A  LINEAR  PROORANMINQ  APPROACH  TO  THE 
CHEMICAL  EQUILIBRIUM  PROBLEM 

1 .  INTRODUCTION 

The  extension  of  linear  progranmlng  methods  into  the  non¬ 
linear  area  Is  recognized  as  one  of  the  outstanding  areas  of 
research  today.  One  way  to  treat  a  general  convex  objective 
function  is  to  locally  linearize  It  by  taking  partial  deriva¬ 
tives.  However  the  authors  believe  that  researchers  In  the 
linear  programming  field  are  Interested  In  Ideas  that  convert 
a  "near”  separable  convex  function  Into  a  completely  separable 
case  where  more  efficient  methods  can  then  taxe  over  (ij,  [2j . 

Although  the  particular  problem  resolved  has  a  non— manage¬ 
ment  application,  it  should  also  be  of  interest  to  a  rather 
large  audience  scheduling  for  tne  Petroleum  Industry.  These 
are  for  the  most  part  chemical  engineers  who  are  familiar 
with  the  application  to  the  Chemical  Equilibrium  Problem  and 
who  can  tane  advantage  of  the  solution  for  their  respective 
companies.  They  will  find  the  companion  paper  [}]  discusses 
more  of  the  chemical  background  to  the  problem  and  gives  an 
alternative  procedure.  The  present  paper  stresses  the 
mathematical  development  of  the  linear  programming  approach 
to  this  problem. 
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The  determination  of  the  chemical  composition  of  a  con4— 
plex  mixture  under  chemical  equilibrium  conditions  is  a 
classic  problem,  l^ere  have  been  many  computing  techniques 
proposed  and  the  constant  appearance  of  new  ones  attests  that 
none  are  entirely  satisfactory.  In  our  second  paper  [j>]  , 
entitled  "Chemical  Equilibrium  in  Complex  Mixtures",  the 

I 

formulation  of  two  methoas  discovered  in  the  course  of  our 
researches  were  reviewed,  one  of  which  is  a  steepest  descent 
based  on  a  quadratic  fit  to  the  free  energy  function  to  be 
minimized,  and  the  other  which  reduces  it  to  a  linear  program¬ 
ming  problem,  purpose  of  this  paper  is  to  give  a 

complete  account  of  the  latter.  Our  purpose  is  to  show  that 
there  is  an  elegant  way  to  transform  the  free  energy  function 
Into  a  convex  separable  function  which  pemits  convenient 
piece— wise  linear  approximation  and  consequent  solution  by 
linear  programming.  Any  desired  degree  of  accuracy  can  be 

reached  by  successively  Improving  the  approximation. 

I 

2.  THE  PROBLEM 

1 

We  consider  an  equilibrium  mixture  containing  m  diff— 

I  erent  atom  types.  While  in  theory  these  will  combine  into 
all  chemically  possible  molecular  species,  in  practice  only 


standard  types  are  considered  including  the  monotonic  types 


p-ioi:o 


which  are  known  to  occur  in  measurable  amounts . 


Let  •  the  number  of  atomic  weights  of 

species  1  present  In  the  mixture. 

Xj  ■  the  nvimber  of  moles  of  molecular  species 
J  present  in  the  mixture  where 

(1)  J  -  1,  2,  ...  n. 

"x  ■  the  total  number  of  m.oles  of  gas  In 
the  mixture,  l.e. 

(2) 

the  number  of  atoms  of  species  1  In 

a  molecule  of  species  J. 

Then  the  mass  balance  equations  are 

n 

(3)  Z  ®1J^J  •  ^1  ^*01'  1  -  1,  2,  . . .  m. 

Th.e  determination  of  thie  equilibrium,  composition  of  a 
gaseous  mixture  Is  equivalent  to  the  determination  of  the 
values  of  the  mole  numbers  Xj  that  obey  constraint  (^)  and 
minimise  the  total  free  enei’gy  of  the  mixture  given  by 


n  n 


which  can  be  shown  to  be  a  convex  function.*  The  values  Cj 


[3]  contains  a  direct  proof;  alternatively  the  fact 
that  the  chemical  equilibrium  problem  can  be  reduced  to  a 
linear  programming  problem  to  any  desired  degree  of  accuracy 
also  proves  convexity. 


/«). 
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are  the  modified  Oibba  free  energy  function  f^/HT  of  the  atonic 
species  at  a  given  temperature  plus  the  natural  logarithm  of 
the  pressure  in  atmospheres. 

Our  problem  is  to  minimise  (4)  subject  to  the  linear 
equality  and  Inequality  constraints  (l),  (2),  (3)* 

In  order  to  apply  linear  programming,  we  make  a  piece¬ 
wise  linear  approximation  to  each  of  the  terms  (xj/x)  In  (Xj /x) 
that  appear  in  (4).  if  we  set  a  -  P  -  a  In  a  then 

we  shall  replace  each  such  curve  by  a  broken  line  function 
0  »  P(a)  such  ao  the  one  below. 


/3 


The  k  points  where  the  two  curves  agree  are  denoted  by 
where  1  »  1,  2,  ...»  k. 

Let  us  imagine  that  the  values  of  x  and  x^  are  fixed  for 
the  moment  and  that  the  values  of  2  0,  x^^g  ^  ®'***'*ln  ^  ® 
are  chosen  so  to  satisfy 


(5) 

(6) 
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Ik 


and  to  minimize  where 

(7)  Zj  -  ♦  ^2^12  ■♦■  •••  "*■  Pk’^ik* 

We  shall  now  prove  that  P  «  P(a)  Is  convex,  that 
(S)  Min  •  X  3(x^/llc)  *  x^  ln(x^/x). 

Proof:  It  Is  clear  that  3  -  p(a)  will  be  convex  If  0 
is  convex  and  this  follows  by  noting  that 

d£ 


alna 


(9) 


■  1  +  In  a 


is  monotonically  increasing.  Next  let  us  substitute 


(10) 

^11  " 

*12 

•  X2X, 

•• •»  *lk  " 

then  (5), 

(6) ,  and  (7)  may  be 

rewritten 

(11) 

1 

X2  + 

X^ 

(12) 

(x^/»^  -  Oj 

[\  ^ 

OjXj  + 

.  . .  +  o^X^ 

(13) 

{z^/f.)  - 

[\ 

SjXg  + 

...  4-  Pj^Xj^ 

and  the  problem  Is  equivalent  to  flnllng  >0,  X^  >  > 

satisfying  (ll),  (12)  for  fixed  x^^  and  x  minimizing  (z^/x). 

If  we  Interpret  X.  >  0  as  the  weights  assigned  to  the  points 
k  ^  “  k 

(aif^i),  then  ^OjX^  and  ^^X^  are  coordinates  of  the  center 
of  gravity  of  the  points.  Hence  we  are  seeking  weights  to 
assign  to  the  points  such  that  the  abscissa  of  the  center  of 
gravity  Is  see  (l2),  and  the  ordinate  ,  see  (13), 

Is  as  small  as  possible.  Obviously  this  smallest  value,  for 
any  convex  curve  ^  ■  P(a\  Is  p(xj/3c),  emd  this  value  Is  ob— 
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talned  for  a  broken  line  function  by  assigning  -  0  to  all 
points  except  the  two  points  on  either  side  of  (x^/x)  and 
weighting  up  these  two  points  appropriately. 

To  solve  the  chemical  equilibrium  problem  by  linear 
programming,  consider  the  problem  of  finding  Xj  >  0,  Xjj  0 
tnd  minimum  z  satisfying 


^11^1  *  *  ^in^n 

®21^1  * 

■  ‘>1 
-bg 

*ml^l  ^mn^'n 

X.  •«■...  -f  X. 

1  n 

—  X 

-  0 

k 

^11 

-  X  X., 

-  0 

-  0 

—  X 

o 

1 

•• 

k 

"  ^1 

«  0 

-  *2 

»  0 

Ic 

“  ^n 

K  K 

K 

"1*1  ♦  "n*n 

•  •  •  ■  * 

It  will  now  be  easy  to  see  that  the  values  Xj  ■  Xj,  that 
form  part  of  the  optimal  solution  to  (14),  give  the  optiaal 
solution  to  (l),  (2),  (3)  If  (4)  Is  replaced  by  the  approxi¬ 


mation 


(15) 
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Proof:  Consider  an  optimal  solution  to  (14)  and  let  the 
values  cf  •  x®.  For  these  fixed  values  of  Xj,  it  Is  clear 
from  the  structure  of  (l4)  that  the  optimal  choice  of  values 
^11*^12*  *'* '*lk  example  must  satisfy  (5)  amd  (6)  and 
minimize  given  by  (7).  Moreover  by  (8)  the  value  of 
Zj  -  x^p(xj/X°)  cmd 

Min  2  +lx°e(x5/x®)  .  V° 


Conversely,  consider  an  optimal  solution  to  (l), 
with  values  Xj  »  xj,we  can  always  find  values  Xj^ 
such  that 


(2),  (3),  (15) 


>  0 


As  we  saw  in  the  discussion  following  (13)>  these  x^,  x^^  are 
a  solution  to  (l4)  whose  value  z  =  z*  is  the  same  as  the 
value  Min  P  obtained  by  substituting  x.  =  In  (35);  hence 
Min  P  =  z*.  Prom  the  relations  Min  P  =  z*  ^  Mln  z  and 
Min  F  ^  =  Mln  z  follows  Mln  P  =  Mln  z  and  our  proof  Is 

conv)lete. 

This  approach  has  several  advantages  over  previous  methods, 
(l)  It  uses  a  standardized  code  of  the  simplex  method 


of  linear  programming. 
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(2)  Only  one  curve  a  log  a  for  0  <  o  <  1^  to  be 

approximated  and  this  can  be  done  as  accurately  as  desired 
since  the  number  of  equations  Is  not  Increased. 

(3)  There  Is  no  necessity  of  deciding  in  advamce  wb.lch 
are  the  so  called  major  components  and  which  are  the  minor 
ones  as  In  the  usual  case  for  other  methods. 

The  evaluation  of  log  can  be  Included  In  the  coding 
so  that  the  successive  approximations  could  be  carried  out 
automatically  to  any  degree  of  accuracy. 

In  order  not  to  tax  the  memory  of  the  computing  machine 
and  also  to  reduce  the  number  of  computations  of 
3^  «  log  a^,  the  following  system  of  "screening"  should  be 
added  to  the  code. 

First  compute  6  “  1  with 

grid  size  .5.  Then  after  the  first  linear  program!  approxi¬ 
mation  has  been  solved,  for  each  tern  (xj/x)  In  (Xj/i?), 
halve  the  grid  size  amd  compute  only  those  new  values  which 
are  adjacent  to  the  current  value  of  a  •  ^j/**  on  the 

second  piecewise  approximation,  if  some  Xj  •  0,  compute 

discarding  the  value  at  ■  1; 

If  Xj/x  -  .5»  compute  log  at  .25  and  .75  and  discaixi  the 
values  at  0  and  1.  If  Xj/l?  la  a  weighted  average  of  two  grid 
points  0  and  .5,  then  Include  a  grid  value  of  at  .25  and 
discard  the  value  at  1,  etc.  In  this  way  for  each  new  piece¬ 
wise  approximation  we  have  at  most  three  values  of  and 
log  for  each  J  ■  1,  2,  ...,  n  such  that  the  range  of 
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values  l8  halved  each  tine.  It  has  been  observed  empirically 
that  successive  values  of  *^ave  stayed  within  the  ranges 

set  up  for  all  previous  cycles  l.e.  the  accuracy  has  doubled 
with  each  successive  approximation. 

After  a  sufficient  number  of  cycles  have  been  carried 
out,  any  Xj  which  Is  still  at  zero  level  In  the  approximate 
solution  could  be  adjusted  to  a  trace  level  by  standard 
methods  used  to  get  the  order  of  n^agnltude  of  such  components 
when  the  non— trace  level  components  have  teen  determined  very 
closely . 

Wyien  there  are  a  lerge  number  of  probletrs  to  be  run,  It 

Is  recommended  that  the  simplex  code  be  modified  to  generate 

# 

Internally  In  the  machine  the  columns  associated  with  the 
variables  this  Is  done  It  Is  probably  best  to  deter¬ 

mine  the  best  choice  of  a^,  analytically  for  each  J 
rather  than  by  the  above  grid  technique.  The  present  RAND 
code  Is  being  modified  along  these  lines. 

There  Is  also  another  Improvement  possible  which  takes 
further  advantage  of  the  structure  by  partitioning  the  basis 
Into  two  parts — the  first  part  associated  with  the  first  m 
rows  and  the  other  associated  with  the  remaining  rows.  Thds 
results  In  an  essential  m  x  m  subbasls  whose  Inverse  Is  needed — 
the  remainder  of  the  computation  Is  tf en  carried  out  Implicitly. 
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